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ABSTRACT 

A  system  of  approximate  two-dimensional  equations  for  the  extensional  and  flex¬ 
ural  motion  of  thermopiezoelectric  plates  subject  to  large  electric  fields  is  derived 
from  the  three-dimensional  equations  of  thermopiezoelectricity  for  small  strains 
small  temperature  variations,  and  strong  electric  fields.  The  two-dimensional  equa¬ 
tions  are  derived  by  introducing  an  appropriate  expansion  for  the  mechanical  dis¬ 
placement,  temperature  field,  and  electric  potential  in  the  thickness-coordinate  and 
integrating  the  balance  laws  and  constitutive  relations. 


I.  INTRODUCTION 

Recently  electroelastic  equations  for  the  extensional  motion  of  very  thin  plates 
with  fully  electroded  major  surfaces  subject  to  large  driving  voltages  and  undergoing 
small  strain  were  derived  [1]  from  the  general  nonlinear  electroelastic  equations  [2,3]. 
These  equations  have  been  used  to  analyze  of  static  and  dynamic  deformations'  of 
laminated  beams  and  plates  due  to  piezoelectric  actuators  [4-6].  A  set  of  two- 
dimensional  equations  was  subsequently  obtained  [7]  for  a  higher  order  approximate 
description  of  the  mechanical  and  nonlinear  electrical  behavior  of  the  relatively  thin 
electroelastic  plates  in  extensional  and  flexural  motion.  The  equations  are  valid  for 
either  the  electroded  or  unelectroded  plates  and  hence  have  been  used  [8]  to  analyze 
partially  electroded  actuators  which  have  been  shown  to  have  a  distributed  shear 
stress  rather  than  a  singular  shear  stress  of  fully  electroded  actuators.  This  prevents 
the  delamination  of  the  piezoelectric  layer  from  the  base  plate. 

In  the  above  analysis,  only  piezoelectric  effect  is  considered,  and  thermal  effect 
is  not  included.  According  to  a  recent  review  article  by  Rao  and  Sunar  [9],  the 
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composite  intelligent  structures  have  different  response  characteristics  at  different 
temperatures  and  the  temperature  variation  in  the  piezoelectric  materials  can  affect 
the  overall  performance  of  the  control  system,  hence  thermal  effects  are  important 
in  the  precision  distributed  sensing  and  control  of  intelligent  structures.  Accord¬ 
ing  to  Rao  and  Sunar ,  the  applications  of  thermopiezoelectricity  theory  to  practical 
engineering  problems  in  general  and  vibration  control  of  flexible  structures  in  par¬ 
ticular,  are  very  few  in  the  literature  and  the  development  of  thermopiezoelectric 
sensors/actuators  is  important  for  advanced  intelligent  structures. 

Mindlin  [10]  derived  a  set  of  two-dimensional  plate  equations  based  on  the  linear 
theory  of  thermopiezoelectricity  for  small  strains,  small  temperature  variations,  and 
weak  electric  fields.  Recent  works  on  laminated  thermopiezoelectric  plates  [11  12]  are 
also  for  the  linear  theory  and  weak  electric  fields.  It  has  been  noted  that  piezoelectric 
materials  are  often  operating  under  large  driving  voltages  or  strong  electric  fields 
[1]  hence  nonlinear  terms  in  electric  fields  should  be  considered.  Besides,  in  [11,12] 
only  equations  of  the  balance  of  linear  momentum  are  derived  but  the  electric  charge 
equation  and  the  heat  equation  are  not  included.  The  electric  and  thermal  fields 
have  to  be  obtained  separately  and  hence  these  equations  can  not  be  applied  to 
structures  with  partially  electroded  actuators. 

A  set  of  equations  for  electroded  very  thin  thermopiezoelectric  plates  subject  to 
large  driving  voltages  or  strong  electric  fields  in  extensional  motion  is  derived  in  [13] 
based  on  the  general  nonlinear  equations  of  thermopiezoelectricity  [14];  the  work  in 
[13]  generalizes  that  reported  in  [1]  to  include  thermal  effects.  Here,  a  system  of 
two-dimensional  equations  for  the  extensional  and  flexural  motion  of  thermopiezo- 
electnc  plates  subject  to  large  electric  fields  is  derived  from  the  three-dimensional 
equations  of  thermopiezoelectricity  for  small  strains,  small  temperature  variations, 
and  strong  electric  fields.  The  two-dimensional  equations  are  derived  by  introducing 
an  appropriate  expansion  for  the  mechanical  displacement,  temperature  field,  and 
electric  potential  in  the  thickness-coordinate  and  integrating  the  balance  laws  and 
constitutive  relations  through  the  thickness.  The  resulting  equations  are  reduced 
to  the  uncoupled  system  of  equations  describing  extensional  motion  and  elementary 
flexure.  These  equations  generalize  those  in  [7]  by  including  thermal  fields.  The 
equations  are  valid  for  either  the  electroded  or  the  unelectroded  plate  and  hence  can 
be  used  to  analyze  either  fully  or  partially  electroded  actuators. 


II.  BASIC  EQUATIONS 

It  has  been  shown  [13]  that  the  thermoelectroelastic  equations  for  infinitesimal 
strains,  smaU  temperature  variations,  and  strong  electric  fields  may  be  written  as 


■^>7,1  =  P^j 
Di,i  =  0 

-  Qi.i  =T pi] 

~  ^ijklSkl  —  ^kijEk  —  Xijd  —  -bkhj  EkEl 
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A-  =  CiktSki  +  €ikEk  +  PiO  +  ^XkjiEkEj 

PV  =  +  pkEk  +ad  +  ^■jkjEkEj 

Qi  —  ^ij^J 

Ski  =  -^{y-kj  +  ut^k) 

A  =  -<l>,k  ■ 

Here  a  comma  followed  by  an  index  /  denotes  partial  differentiation  with  respect 
to  the  referential  coordinate  x;,  a  dot  over  a  variable  denotes  partial  differentiation 
with  respect  to  time  and  repeated  indices  are  to  be  summed.  The  range  of  indices 
i,j,  k,  l  IS  1, 2, 3.  In  Eq.  (1)  p  is  the  mass  density,  r.y  the  stress  tensor,  u.  the  mechan¬ 
ical  displacement,  Sij  the  infinitesimal  strain  tensor,  Z),  the  electric  displacement 
vector,  <f>  the  electric  potential,  Ei  the  electric  field,  Qi  the  heat  flux,  t?  the  entropy 
density,  T  a  uniform  reference  temperature,  and  $  the  temperature  variation,  c.yi-;, 
e.t/,  denote  the  elastic,  piezoelectric  and  dielectric  constants,  A,_,  the  thermal 
elastic  constants,  p,  the  pyroelectric  constants,  a  is  related  to  specific  heat  and  «:,y 
is  the  heat  conductivity  tensor.  S*/,y,  Xkji,  and  y.y  are  nonlinear  material  constants. 
We  note  that  ‘bkuj,  the  effective  electrostrictive  constants,  include  the  effect  of  the 
Maxwell  electrostatic  stress  tensor  [13].  Equation  (1),  is  the  balance  of  linear  mo¬ 
mentum,  (1)2  the  electrostatic  charge  equation,  (1)3  the  heat  equation,  (1)4_7  the 
constitutive  equations,  (1)3  the  strain-displacement  relation,  and  (1)9  the  electric 
field-potential  relation.  When  the  electric  fields  are  weak,  quadratic  terms  on  the 
right-hand  side  of  equations  (1)4.6  are  dropped  and  Eq.  (1)  then  reduces  to  that  for 
linear  thermopiezoelectricity. 


III.  TWO-DIMENSIONAL  EQUATIONS  FOR  THIN  PLATES 

A  plan  view  and  a  cross-section  of  the  thin  plate  with  thickness  2h  are  shown 
in  Figure  1  along  with  the  coordinate  system.  The  top  and  bottom  surfaces  may 
be  identically  electroded;  the  electrodes  in  a  given  region  are  assumed  to  have  a 
potential  difference  or  voltage  across  them. 

We  now  obtain  approximate  two-dimensional  mechanical  equations  of  motion, 
charge  equation  of  electrostatics  and  heat  equation  for  the  thin  plate  shown  in  Figure 
1  by  employing  assumed  expansion  of  uy,  and  0  in  the  thickness  coordinate. 
Since  only  the  lowest  frequency  approximations,  i.e.,  equations  of  the  extensional 
and  flexural  motion  of  plates  are  of  interest  here,  we  expand  uj  in  the  form 

Uj  =  uf\xa,t)  -f  X3uj'*(Xa,  t)  +  xluf\xa,  t)  (2) 

in  which  Uj  is  included  to  allow  for  the  free  thickness-strains  accompanying  ele¬ 
mentary  flexure  [7].  We  also  adopt  the  convention  that  indices  a,6,c,d  take  values 
1,2.  Considering  the  behavior  of  ^  across  the  electrodes  when  they  are  present, 
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Fig.  1,  Plan  view  and  cross-section  of  a  partially  electroded  thin  plate. 


we  expand  <f>  in  the  form 

^  +  (^  -  +  f  (fl  -  (3) 

in  whia  only  the  4,^  tettn  contribute  to  the  volt.ge  across  the  electrodes.  Equation 
(3)  basically  taken  from  [7],  with  the  difference  that  here  we  have  but  the 
corresponding  equation  in  [7]  does  not;  ^(<>)  is  needed  for  a  complete  description  of  </, 
n  the  une  ectroded  region,  although  in  the  case  treated  in  [8]  it  has  no  contribution^ 
In  an  electroded  region  the  two-dimensional  plate  electric  potentials  ^«»  and  6^ 
must  be  independent  of  coordinates  in  the  plane  of  the  plate.  However  in  an 
unelectroded  region  and  <l>(  )  are  in  general  functions  of  coordinates  x„.  Clearly 

ceded  to  fuUy  describe  the  response  of  a  partially  electroded  thin  plate.  For  the 
temperature  field  6,  we  take  [10]  ^  ® 

0  =  d^°'>{Xa,t)  +  Xie^^'>(Xa,t)  (4) 

bv  tTm  TAt’";  TTAr  ‘Charge  equation  (l)^ 

y  {  ,X3/2h,x3/h  -l,ixl/h^-l)x3/h},  the  heat  equation  (1)3  by  {l.xa},  integral 

.he  result, „g  „u..„„s  from  -4  to  i  .„d  obt.iu  the  lollowing  „i„e  teo-dimenSul' 
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equations  corresponding  to  the  balance  of  linear  momentum 


rSa  -  nr^r’  +  Fr  =  P  E  ,  n  =  0. 1, 2 


J«-i) 


?(") 


m=0 


four  two-dimensional  electric  charge  equations 

j^i>a-^DS"'+4'>=o 


i-£){2)_  n(0)_  An(i)_n 

-L/jW  _  in(i)  _  An(2)  ,  in(o)  _  n 


/,3-a.a  ^3 

and  two  two-dimensional  heat  equations 


(5) 


(6) 


«?]  +  4“>  = 

<?i|l  -  Q?  +  4'*  =  -rpi7<'* 


where 

i'^hi  ^  xjlTki,  Dk,Qk,  'n]dxz 

F'r^  =  KrayjE  ,  4"^  =  ■  4"’  =  [xJQalE 

i/00  =  2h  ,  Hqi  =  iJio  =  0  ,  =  \h^ 

O 

Ho2  =  H20  =  -h^  ,  Hi2  =  H21  =  0  ,  H22  = 

*>  5 


(8) 


and  we  note  that  terms  Fj^\  ^3”^  and  ^3”^  arise  from  the  integration  by  parts. 

Terms  ^3  ^  and  ^3  \  which  are  analogous  to  ^3*^^  and  vanish  because  the  X3 
dependence  of  coefficients  of  and  vanish  at  X3  =  ±h.  We  also  note  that 
in  the  balance  of  linear  momentum  equations  (5)  the  three  equations  for  n  =  2 
will  not  actually  occur  because  they  will  be  eliminated  by  allowing  for  the  free 
development  of  accompanying  anisotropic  flexure  [7].  In  an  elect roded  region 
in  which  ~  V  where  V  equals  the  voltage,  is  a  constant  and  the  first  two 
equations  of  (6)  are  not  needed  to  obtain  the  solution.  They  simply  serve  to  define 

^3  ^  and  ^3^^  in  terms  of  the  electric  displacement  D.  However,  in  an  unelectro ded 
region  in  which  V  is  not  prescribed  these  equations  are  required  to  obtain  a  solution. 
Similarly,  if  temperature  is  known  on  both  major  surfaces  of  the  plate,  and 
are  determined  directly  from  (4)  and  equation  (7)  is  not  needed  to  obtain  the 
solution;  it  serves  to  determine  93^^  and  in  terms  of  other  variables.  However,  if 
heat  flux  is  prescribed  at  either  one  of  the  major  surfaces  of  the  plate,  equation  (7) 
is  needed  to  obtain  a  solution. 
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At  this  point  we  have  nine  two-dimensional  stress  equations  in  (5).  However, 
since  we  are  interested  in  obtaining  only  the  uncoupled  equations  of  anisotropic 
extension  and  elementary  flexure,  it  is  convenient  to  rewrite  these  as  separate  es¬ 
sentially  extensional  and  essentially  flexural  equations.  We  note  these  equations  are 
coupled  due  to  anisotropy  of  the  plate  material.  We  further  note  that  we  can  not 
complete  the  reduction  without  the  plate  constitutive  relations,  which  are  obtained 
m  the  next  section.  For  the  plate  constitutive  equations  it  is  convenient  to  define 
Mmdlin’s  plate  strains,  which  is  done  simply  by  substituting  the  expansion  (2)  into 
the  strain-displacement  relation,  Eq.  (1)3,  and  rearranging  terms;  the  result  is 


n=:0 


4’  =  jK‘? + -S’ + (» +  +  {.,-!"«>)] 


(«) 


(9) 


We  now  note  that  we  will  not  need  s\f  because  of  the  reduction  that  is  to  be  made. 
When  elementary  flexure  and  extension  are  to  be  uncoupled,  from  (5),  following  the 
steps  in  (7],  we  may  separate  the  essentially  extensional  plate  equations. 


(10) 


.(2)  . 
^ab,a 


^3i 


r(2) 


from  the  essentially  flexural  plate  equations. 


2 

3' 

^ab,a  ^3b  +^b  =  ^P^  H 


(11) 


I  -t- 

‘aZ,a  -^^33  '  ^3 


the  case  of  extension  we  must  first  allow  for  the  free  plate  thickness-strains 
^33  by  setting  r33  =  0.  Then  in  order  to  eliminate  the  first  order  extensional 
equations  completely  we  set  =  0.  This  has  the  effect  of  eliminating  the  second 
order  equations  in  (10)  completely  since  all  second  order  plate  stress  resultants 
are  irrelevant  in  the  approximation,  and  may  be  ignored,  as  may  the  dynamic  ter^^is 
on  the  right-hand  side  of  (10)3  and  (11)3.  Furthermore,  in  order  to  eliminate  flexure 
from  the  extensional  equations,  from  (11)  it  is  clear  that  we  must  have  =  0 
=  0,  7-33  =  0.  Collecting  all  the  conditions  on  the  stress  resultants,  we  have 


M  - 


=  0 


=  0 


(12) 


the  first  of  which  will  be  used,  in  the  next  section,  to  reduce  the  general  plate 
thermopiezoelectric  constitutive  equations  to  those  that  are  suitable  for  anisotropic 


660 


extension.  Since  we  have  eliminated  Hexure,  therefore  =  0.  Also  since  we  are 
well  below  the  lowest  thickness  resonant  frequency  of  the  plate,  we  may  assume 
Uj  =  0,  iuj  =  0,  then  all  that  remains  out  of  (10)  is 


=  2phul 

which  are  the  two  equations  for  the  extensional  motion  of  thin  plates. 

(^)For  flexure,  we  must  first  allow  for  the  plate  thickness  strain  by  setting 

r33  =  0.  This  has  the  additional  effect  of  eliminating  the  second  order  flexural 
equation  completely  because,  as  already  noted,  all  second  order  plate  stress  resul¬ 
tants  are  irrelevant  in  this  approximation,  and  may  be  ignored.  Then  in  order  to 
eliminate  extension  from  the  flexural  equations,  from  (10)  it  is  clear  that  we  must 

have  =  0,  =  0,  =  0.  These  conditions  enable  us  to  write 


.•■.(0) 


(13) 


.(«)  - 


=  0 


JO) 


=  0 


(14) 


the  first  of  which  will  be  used,  in  the  next  section,  to  reduce  the  general  plate 
thermopiezoelectric  constitutive  equations  to  those  that  are  suitable  for  elementary 
flexure  of  thin  plates.  Since  we  have  eliminated  extension,  therefore  =  0.  Also, 
since  we  are  well  below  the  lowest  thickness  resonant  frequency  of  the  plate,  we  may 
assume  that  tiy  =  0,  uf'  =  0.  In  order  to  complete  the  reduction  to  the  elementary 

theory  of  flexure  we  must  take  the  thickness  shear  strains  51°^  to  vanish  17]  with 
which  (9)2  yields 

(15) 


which  enables  us  to  obtain  a  single  equation  in  the  one  dependent  variable  in  the 
elementary  theory  of  the  flexure  of  thin  plates.  Utilizing  ^  q,  which  eliminates 
rotary  inertia,  we  obtain  rjj  from  (11)2,  which  when  substituted  into  (ll)i  yields 

^lb!a6  +  ^b!b  +  ^3°^  =  (16) 

This  is  the  equation  for  the  elementary  theory  of  the  flexure  motion  of  thin  plates. 

Thus,  at  this  stage  we  have  the  extensional  equations  of  motion  (13)  which 
require  constitutive  equations  for  the  flexural  equations  of  motion  (16)  which 
need  constitutive  equations  for  the  four  plate  equations  (6)  for  electrostatics 
which  necessitate  constitutive  equations  for  (n  =  0, 1,2,3),  and  two  plate  heat- 
conduction  equations  (7)  which  require  constitutive  relations  for  <5^"^  (n  =  0, 1)  and 
i/l")  (n  =  0, 1).  As  noted  earlier,  the  constitutive  equations  will  be  derived  in  the 
next  section  in  terms  of  the  nine  dependent  variables  </J")  (n  =  0, 1  2  3)  and 
^(")(n  =  0,l).  ' 

To  the  foregoing  plate  equations  we  must  adjoin  the  appropriate  initial  and 
boundary  conditions.  For  the  equations  of  anisotropic  extension  (13),  at  an  interface 
separating  one  region  from  another  we  have  the  well-known  continuity  condition 

'^alri“^j  =  0,  lui'’’]=0  (17) 
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where  ria  denotes  the  unit  normal  directed  from  the  to  the  “-f”  side  of  the 
interface  and,  in  which  we  have  used  the  usual  notation  [Aj  ^  for  the  jump. 

At  the  edge  either  or  or  some  combination  thereof  is  prescribed.  For  the 
equation  of  elementary  flexure,  at  an  interface  we  have  the  continuity  condition 


ds 


where 


^n3  —  ^6^63 


1  =  0, 

1^1=0,  [.(“'1=0 

(18) 

tW 

‘ns 

(1) 

(19) 

and  sij  denotes  a  unit  vector  tangent  to  an  interface  of  separation  in  the  counterclock 
direction.  At  an  edge  either  and  T^^^+dri^J /ds  or  dJ°^/dn  and  or  some 

combination  thereof  is  prescribed.  In  addition,  there  are  the  well-known  conditions 
across  corners  of  discontinuous  curves. 

In  the  case  of  the  two-dimensional  plate  equations  of  electrostatics  (6),  at  an 
interface  separating  one  region  from  another,  we  require  the  continuity  conditions 


[^Wj  =  0,n  =  0,l,2,3  (20) 

In  addition,  by  multiplying  the  well-known  three-dimensional  electrostatic  continuity 
conditions  10^7x4  =  0  by  {l,x^/2h,xl/h^  -  l^xj/k^  ^  l)x3/h},  and  integrating 
the  resulting  equations  from  -h  to  A,  we  obtain  the  following  two-dimensional  plate 
continuity  conditions 


nalD^°^j  =  0  ,  =  0 

=  0  ,  n4Z)W-h2Z)(i)]  =  0 

At  an  edge  either  -  K^D^a^)  and  or 

<jA^\  and  (fP)  or  some  combination  thereof  is  prescribed. 

For  the  heat  equations  (7),  the  continuity  equations  are  [10] 

[^'^“^1  =  0  ,  =  0  ,  n„(g<“)|  =  0  ,  n4gW]  =  0  (22) 

At  an  edge  either  6^^^  and  or  UaQ^a'^  and  riaQ^a^  or  some  combination  thereof  is 
prescribed. 


IV.  TWO-DIMENSIONAL  CONSTITUTIVE  EQUATIONS 

In  this  section  we  obtain  constitutive  relations  for  the  stress  tensor,  electric  dis¬ 
placement,  heat  flux  and  entropy  defined  in  (8).  The  resulting  constitutive  equations 
are  then  reduced  to  those  appropriate  for  the  uncoupled  equations  of  anisotropic  ex¬ 
tension  and  elementary  flexure  by  employing  Eqs.  (12)i  and  (14)i,  respectively. 
Clearly,  the  constitutive  equations  are  obtained  by  substituting  from  (1)4-7  into 


662 


(8)i,  employing  (9)  and  the  analogous  equations  for  electrical  equations  and  in- 
tegraUng  through  the  thickness.  From  (1)9  and  (3)  we  find  that  the  analogous 
electrical  equations  take  the  form 

As  noted  earlier,  in  this  low  order  treatment,  we  ignore  the  second  order  equa¬ 
tions  completely,  and  we  do  not  need  the  plate  strains  in  (9)i,  which  means 
that  for  our  purposes  here,  Eq.  (9)i  takes  the  form 

^ij  =  ^  (24) 

T»=0 


which  is  analogous  to  (23). 

For  convenience,  we  introduce  the  usual  compressed  matrix  notation  for  stresses 
and  strains.  In  this  convention  the  tensor  indices  ij  or  kl  are  replaced  by  p  or  ^ 
which  take  values  1,  2,  3,  4,  5,  6  as  ej  or  kl  take  values  11,  22,  33,  23  or  32,  31  or 
13,  12  or  21  respectively.  Accordingly,  we  write  constitutive  equations  (1)4-6  in  the 
compressed  notation  as 


Tp  —  CpqSq  —  ekpEk  —  Xp6  —  -bkipEkEi 

Di  =  CiqSq  +  tikEk  +  PiO  +  ^XkjiEkEj  (25) 

PP  =  XqSq-{-  pkEk  -f  -f-  -'JkjEkEj 

where  the  sum  from  1  to  6  on  the  repeated  matrix  indices  is  understood.  We  also 
note  that  in  the  matrix  notation  Eq.  (24)  takes  the  form 

=  X]  (26) 

n=0 


Similarly,  from  Eq.  (8)1  we  obtain 

^  ^  J  ^  ^  I  h 

Now,  substituting  from  (25)  into  (27)  and  recalling  the  definitions  of  and 

T]^ ^  in  (8)1,  substituting  from  (26),  (23),  (4)  and  (1)7,  and  integrating  through  the 
thickness,  we  obtain  the  following  plate  constitutive  relations 

rW  =  2Acp,5f)  +  -  2AA,0(°)  - 


Df^  =  2heu,Sf)  -  2he,^<i,f  +  1a£,,«5«  -  £*3^(1)  +  2hpte(^)  +  ^x.y^N^ 

4'^  =  -fi^eipSp^  -  ^h^€ta<^(p  +  Aa2£^__^(3)  _  1a£43^(2)  +  +  ^Xijt^jp 

=  ^A^e,p5W  -  +  ±h^e,^^(p  -  i/.2£,3^(i) 

-  +  ixoiAfJ)  (28) 

^  =  5^®«*p'S'F^  -  ^h^^ka<f>P  +  ^h*ekapp  -  ^A^£i3fli(2)  + 

/>r;«»  =  2/tAp5(°)  -  2/.p<.^{»)  +  hp^0  -  p3^(i)  +  2A£.^«»  +  i7.,7v/«) 

=  ^h^^pSp'>  -  -h^p^4,p  +  ^h'^p^^ip  _  ^hp3<f>^^'>  +  ^h^a0^^'>  + 

Qp^  =  -2hK,j{$f  +  ^3y0(l)) 

== 

where  the  nonlinear  electrical  terms  N^p  (n  =  0, 1,2,3)  are  defined  by 

Kl'  =  '■P^SV?’  +  ^4”/."  -  j(4’4”’ + 4’4">) 

-  s(4‘'4’ + 4”4")  + 1|4>4'’  +  54>4>i 
=  (^(*'”4’  -  '*'“’4’)  + 1(4>4'  -  ■#"’4’)  +  *<‘'4’i 


45 

4? = 4 


4?  =  ‘’in|(4'4> + 4>4')  -  5(4”4’ + 4’4"') 

- 1(4>4’ + 4’4‘') + i(4>4> + 4)4')] 

4;>  =  Mi^'«4>  +  44>4>  -  4^«4'  -  Uyp  -  4^(.)^P)  + 


4?  =  (i4'*<‘>  +  H*p>4>i 

»r('2)  ,  If  J-D  .(2)  4  ,  //■ 


4?  =  ‘’fe4’4’  -  ^(4’’'?’ + 4’4')  - 1(4’4'  +  4'4')  ‘  ’ 

^iS4>4’  +  44)4'  +  |4)4'i 

4?  =  »’i5|4'4>  +  g«p>4'  -  ^^04)  +  i^w4>  - -  j|*(»)4>| 
4?  = 
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= ‘‘I  + /;vg') 

-  |wg'/.'’ + +  /."'^g>)l 

Since  we  are  interested  in  the  low  frequency  range  of  extensional  and  flexural 
motion,  we  do  not  need  the  shear  and  thermal  correction  factors  [10]  in  the  con¬ 
stitutive  relations.  These  constitutive  relations  are  not  yet  in  a  useful  form  for  our 
purposes  because  we  have  not  yet  imposed  the  extensional  plate  relaxation  condi¬ 
tions  (12)i  and  flexural  plate  relaxation  conditions  which  in  the  compressed 

matrix  notation  take  the  forms 

=  Tf*=0,  (30) 

.i"  =  0,  ,«=0,  r«  =  0  ,3„ 

In  view  of  (30)  and  (31)  we  introduce  a  matrix  index  convention  which  will  be 
of  considerable  use  in  the  sequel.  We  let  subscripts  u,  v,  w  take  the  values  3,  4,  5 
while  subscripts  t  take  the  remaining  values  1,  2,  6.  This  permits  us  to  write 
(30)  and  (31)  as 

=  0  ,  =  0  (32) 

We  may  now  write  (28)i,2  in  the  form 

r(«)  =  2Ac„5(“)  -h  2hc„5(‘»  +  2/ic„^(»)  - 

o 

-t-  e3r<l>^^'>  -  2hX,0(’>)  -  ^hijrNlf 
r(«)  =  2Ac,,,5(“)  +  2Ac^„S^)  +  2Ae„„,^W  -  ^ 

4  9  1 

+  3 

Equations  (33)2,4  rnay  be  readily  solved  for  5"^^^  and  respectively,  with  the 


result 
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^  eau»Cu,*(^W  ^  -eawCy,l4>^l^  -  ^e3u,Cu,^(^(’) 

+  A„c-i^(°)  +  ^b.,^c-lNlp 

_  Ac-^e3^^(2)  ^  ^ 

+  A„c-i^(i)  + 

where  the  matrix  suras  are  over  the  indices  3,  4,  5  as  a  result  of  the  convention 
bubstitution  from  (34)  into  (33)i^3,  yields  ^  ’ 

r/o)  =  2Ar„5i‘>)  +  +  ^3,^(1)  _  2AAP5W  -  IS?- 

’  3  ’  ^  ^  2 

where 

Trj  =  Cr,  -  CrvC^^Ctos  ,  V’ia  =  Ck,  -  ekv,C~lcvs 

A?  =  Aj  -  Au,c„ic„a  ,  6fy^  =  6,ya  -  bijy,c~lc„, 

We  note  that  the  Trj  are  the  Voigt’s  anisotropic  elastic  plate  constants  and  the  ^ 
are  the  effective  plate  electrostrictive  constants.  The  constitutive  relations  (35)  are 
in  the  form  required  for  the  equations  for  anisotropic  extension  (13)  and  elementary 
nexure  (lb).  At  this  point  we  note  that  the  substitution  of  (15)  into  (9)2  yields 

=  -“Si  (37) 

which  IS  an  important  relation  for  the  flexural  equation.  Writing  5^"^  in  the  consti¬ 
tutive  equations  (28)3-8  as  two  separate  terms,  one  containing  and  the  other 

containing  5„"  in  accordance  with  our  values  for  the  indices,  and  substituting  from 
(o4j  into  the  appropriate  equations,  we  obtain 

4"’  =  2**.SW  -  2*a.0«  +  -  C„#»  + 

“  “5“*.^!'’  -  3-C*.4'’  +  -  yfts/*'  + 

4”  =  - 1* V4"' + 

-  +  ji'piS'*’  +  yAi,,4”*  +  ifcrtJvW 

4”  = 
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=  2/iA?5<«)  -  2h-pA(^  +  -  p3^(i)  +  2hae(^)  + 

(38) 


where 

C^>  —  ^kj  +  ^kv^vw^jtv  ?  Ckc  ~  3ejtc  +  hejtvC^^^Cctf; 

C  =  +  7e^vC|;tiecio  ,  x^-j^  =  xijk  +  A,yjt 

A  ?  -1  (39) 

^xjk  -  OijwC^v^kv  ,  Pk=Pk-\~  ^kvCy^K 

cc  —  t  70  “  70 

The  constitutive  equations  in  (38)  are  in  the  form  required  for  use  in  the  plate 
electrostatic  equation  (6)  and  heat  equation  (7).  However,  in  order  to  use  the 
constitutive  equations  (35)  in  the  balance  equations  (13)  and  (16),  they  must  be 
converted  from  matrix  form  back  into  tensor  form.  This  is  accomplished  simply  bv 
replacing  r  by  ab  and  s  by  cd  wherever  they  occur  in  (35).  For  example 


ril'>  =2hr,,,,si°J  +  -  2hXl,9i<^)  - 

=-^^abcdSil^  +  Y^cab<l>^P  -  '^'^cab^P  + 


Equations  obtained  by  substituting  from  constitutive  equations  (40)  and  (38) 
into  (6),  (7),  (13)  and  (16),  along  with  (9)2  for  and  (37)  are  nine  equations  in  the 

nine  dependent  variables  uf\  <!>(’>'>,  (f>W,  <^(2),  ffW  and  0(^),  which  are  nonlinear 

in  the  plate  potentials  By  proper  substitutions  the  boundary  conditions  can 
also  be  written  in  terms  of  the  same  nine  variables  in  each  region. 


V.  SUMMARY  OF  EQUATIONS 
Balance  laws 


rZ  +  Fp^  =  2phuP 
^Sab  +  FS  +  f'^^=2phuP 
^Dpl+dP  =  0 

-L £)(i) L dW  1  j(i)  _  Q 

-L£)(2)_£)(0)_  Ar,(i)_o 

/j2  “.o  ^a.a  /j2  ^3  -  " 


(41) 


667 


(42) 


_Z)(3)  _  Id(\)  _  An(2)  ,  1  n(0)  _  n 

f^L>a,a  ^3^3  +  ^^3  -0 

<3i?i  +  9i"^  =  -Tp,)(®) 

Strain-displacement  relations 

''S’ =  )(«»+«!”') 

sir = -a 

Constitutive  equations 

r„V  =  2Ar„6<;rf5^2^  +  -  2h\’’.0^°^  -  ii?.  . 

^aV  =  +  Y^cab^W  -  ^ 

Z)f ^  =  'ih^kcdS''^  -  2hCtc<l>f  +  +  2ApiM<>)  + 

4'^  =  ~  jCkc4p  +  ~Ckc<i>^^  -  f  a34i(")  + 

=  ~i>kcdSf}  -  ‘^h\kc4>f  + 

-  — )  +  1x0*4'^ 


jit,  3 

PJ7^“^  =  2AA^j5f°^  -  2Apa^(®)  +  ^Ap„^_(|)  -  ^3^(1)  +  2hae^°'>  + 

P7](^^  =  -  |Ap3«i(2)  +  ^/,3^^(i)  ^  i;y..^.yv^i) 

gr>  =  -2A«.,(0y)  +  A3^#)) 
gp)  =  -h^^i/p 

Continuity  conditions 

«44“^I  =  o,  Iup)]  =  o 


(43) 


(■)..._«  1^1  =  0,  l4'’'|  =  o 


"aKi  Inj  =  0  ,  (r^“)  +  — 

I?i<")l  =  0,n  =  0,l,2,3 
«api“)l=0,  na[£>p)j  =  0 


(44) 
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nalZ)(2)  _  /,2z?(0)j  =  0  ,  n4£>W  -  =  0 

[0«»I  =  O,  I^Wl  =  0.  n4g(<»]=0,  n,[Q(»J  =  0 

[0(“>]  =  O,  [^(’)l  =  0,  naigl“)l=0,  na[gW]  =  0 

Definitions  of  material  constants 

^rv^vuj^tys  j  0Jbs  —  ^ks  ^kw^wv^vs 
^w^yjv^vs  )  —  ^ijs  ^ijw^xvv^vs 

Ckj  ^kj  H”  ^kv^vw^jw  J  Ckc  ^  *^^kc  “f"  ^^kv^vw^cw 

C  *”  ^^kc  +  '^^kv^vw^cw  5  X^jk  —  Xijk  +  ^ijk 
^ijk  ^ijw^u)v^kv  j  Pk  Pk  “I"  ^kv^vw^tv 
^  ^w^^fjoy^v  j  'Ti;  ~  7»J  ^ijw^^wv^v 

Ranges  of  indices 


(45) 


ij,k,l :  1,2,3;  a,6,c,d:  1,2;  p, 9  :  1, 2, 3, 4, 5, 6;  r,s,t  :  1,2,6;  u,v,n;:  3,4,5 


(46) 


VL  CONCLUSIONS 

A  set  of  two-dimensional  plate  equations  has  been  derived  for  the  extension  and 
flexure  of  a  thermopiezoelectric  plate  under  strong  electric  fields.  The  equations 
are  nonlinear  in  the  plate  electric  potentials  and  can  be  applied  to  fully  or  partially 
electroded  actuators  or  sensors.  These  equations  include  the  results  of  [1],  [7]  and 
[13]  as  special  cases. 
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